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Abstract 

We  present  the  modeling  and  optimization  of  a  new  hybrid  solar  thermoelectric  (HSTE)  system  which  uses  a  thermosyphon  to  pas¬ 
sively  transfer  heat  to  a  bottoming  cycle  for  various  applications.  A  parabolic  trough  mirror  concentrates  solar  energy  onto  a  selective 
surface  coated  thermoelectric  to  produce  electrical  power.  Meanwhile,  a  thermosyphon  adjacent  to  the  back  side  of  the  thermoelectric 
maintains  the  temperature  of  the  cold  junction  and  carries  the  remaining  thermal  energy  to  a  bottoming  cycle.  Bismuth  telluride,  lead 
telluride,  and  silicon  germanium  thermoelectrics  were  studied  with  copper-water,  stainless  steel-mercury,  and  nickel-liquid  potassium 
thermosyphon-working  fluid  combinations.  An  energy-based  model  of  the  HSTE  system  with  a  thermal  resistance  network  was  devel¬ 
oped  to  determine  overall  performance.  In  addition,  the  HSTE  system  efficiency  was  investigated  for  temperatures  of  300-1200  K,  solar 
concentrations  of  1-100  suns,  and  different  thermosyphon  and  thermoelectric  materials  with  a  geometry  resembling  an  evacuated  tube 
solar  collector.  Optimizations  of  the  HSTE  show  ideal  system  efficiencies  as  high  as  52.6%  can  be  achieved  at  solar  concentrations  of  100 
suns  and  bottoming  cycle  temperatures  of  776  K.  For  solar  concentrations  less  than  4  suns,  systems  with  thermosyphon  wall  thermal 
conductivities  as  low  as  1 .2  W/mK  have  comparable  efficiencies  to  that  of  high  conductivity  material  thermosyphons,  i.  e.  copper,  which 
suggests  that  lower  cost  materials  including  glass  can  be  used.  This  work  provides  guidelines  for  the  design,  as  well  as  the  optimization 
and  selection  of  thermoelectric  and  thermosyphon  components  for  future  high  performance  HSTE  systems. 

©  2011  Elsevier  Ltd.  All  rights  reserved. 
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1.  Introduction 

Efficient  renewable  energy  sources  are  in  significant 
demand  to  replace  diminishing  and  environmentally  harm¬ 
ful  fossil  fuels.  The  combination  of  commercial  and  resi¬ 
dential  buildings  as  well  as  the  industrial  sector  currently 
consumes  72%  of  the  total  energy  in  the  US  (Davidson, 
2005;  Kelso,  2009).  A  significant  portion  of  this  energy  is 
used  in  the  form  of  heat,  such  as  for  hot  water  heaters  in 
homes  and  food  processing  in  industrial  process  heat 
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(IPH)  applications,  where  temperatures  range  from  50  °C 
to  1000  °C  (Desideri  et  al.,  2009;  Demeter  et  al.,  1991; 
Schnitzer  et  al.,  2007;  Lecuona  et  al.,  2009).  The  abundance 
of  solar  energy  promises  efficient  methods  to  meet  the  cur¬ 
rent  heating  and  electricity  needs.  For  example,  hot  water 
heaters  have  been  commonly  used  in  households  (Han 
et  al.,  2010;  Zhiqiang,  2005),  and  large  scale  solar  thermal 
plants  deliver  distributed  electrical  power  to  cities  (Manuel 
et  al.,  2002;  Pacheco,  2001).  In  most  cases,  however,  these 
systems  are  limited  to  providing  either  heat  or  electricity. 
More  recently,  hybrid  systems  with  either  photovoltaics 
or  thermoelectrics  to  generate  both  electrical  power  and 
heat  have  been  investigated.  Hybrid  photovoltaic-thermal 
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Nomenclature 

A 

area  (m2) 

Tc 

condenser  outer  wall/bottoming  cycle  tempera¬ 

C 

solar  concentration  factor  (-) 

ture  (K) 

cp 

specific  heat  capacity  ( J/kg  K) 

Ts,E 

thermosyphon  evaporator  outer  wall  tempera¬ 

Cs 

pool  boiling  constant  (-) 

ture,  thermoelectric  cold  side  temperature  (K) 

g 

gravitational  constant  (N/kg) 

Psat,E 

thermosyphon  evaporator  saturation  tempera¬ 

G 

incident  radiative  heat  flux  (W/m2) 

ture  (K) 

h 

heat  transfer  coefficient  (W/m2  K) 

P sat,  C 

thermosyphon  condenser  saturation  tempera¬ 

hfg 

latent  heat  of  vaporization  (J/kg) 

ture  (K) 

h/R 

evaporator  thin  film  evaporation  heat  transfer 

TSs 

selective  surface/thermoelectric  hot  side  temper¬ 

coefficient  (W/m2  K) 

ature  (K) 

hp.R 

evaporator  liquid  metal  pool  boiling  heat  trans¬ 

Pvl,E 

thermosyphon  evaporator  liquid  vapor  interface 

fer  coefficient  (W/m 2  K) 

temperature  (K) 

k 

thermal  conductivity  (W/m  K) 

Pvl,C 

thermosyphon  condenser  liquid  vapor  interface 

La 

adiabatic  section  length  (m) 

temperature  (K) 

Lc 

condenser  section  length  (m) 

Tw,e 

thermosyphon  evaporator  inner  wall  tempera¬ 

Le 

evaporator  section  length  (m) 

ture  (K) 

Lf 

evaporator  thin  film  length  (m) 

Tw,c 

thermosyphon  condenser  inner  wall  temperature 

L0 

total  thermosyphon  length  (m) 

(K) 

LP 

liquid  metal  pool  height  (m) 

Too 

environment  temperature  (K) 

Lte 

thermoelectric  leg  length  (m) 

x,y 

coordinates  (m) 

M 

molecular  mass  (kg/kmol) 

m 

mass  flux  (kg/m2  s) 

Greek  symbols 

P 

pressure  (Pa) 

a 

absorbtivity  (-) 

Pc 

critical  pressure  (Pa) 

E 

emissivity  (-) 

p, 

liquid  pressure  (Pa) 

t] 

efficiency  (-) 

Pr 

reduced  pressure  (  ) 

tj* 

efficiency  as  kw  — >  oo  (-) 

Pte 

thermoelectric  power  (W) 

Pr 

reflectivity  (-  ) 

<1 

heat  flux  (W/m2) 

P 

density  (kg/m3) 

Qtoss 

emissive  loss  (W) 

<*B 

Stefan  Boltzmann  constant  (W/m  K4) 

Qout 

bottoming  cycle  heat  transfer  (W) 

a 

surface  tension  (N/m) 

Qso/ar 

solar  heat  input  (W) 

a' 

condensation  coefficient  (-) 

rt 

thermosyphon  inner  radius  (m) 

% 

transmissivity  (-) 

ra 

thermoelectric  inner  radius,  thermosyphon  out¬ 

P 

dynamic  viscosity  (Pa  s) 

er  radius  (m) 

5 

liquid  film  thickness  (m) 

rte 

thermoelectric  outer  radius  (m) 

e 

thermosyphon  inclination  angle  from  the  verti¬ 

Ri 

evaporator  wall  radial  thermal  resistance  (K/W) 

cal  (°) 

R2 

evaporator  liquid  pool  thermal  resistance  (K/W) 

r3 

evaporator  thin  film  thermal  resistance  (K/W) 

Superscript 

r4 

vapor  core  axial  thermal  resistance  (K/W) 

DC 

direct  current 

Rs 

condenser  film  thermal  resistance  (K/W) 

HSTE 

hybrid  solar  thermoelectric  system 

r6 

condenser  wall  radial  thermal  resistance  (K/W) 

IPH 

industrial  process  heat 

Ri 

condenser  wall  axial  thermal  resistance  (K/W) 

PV 

photovoltaic 

R& 

liquid  metal  thin  film  evaporation  interfacial 

PVT 

hybrid  photovoltaic  thermal 

thermal  resistance  (K/W) 

ss 

selective  surface 

r9 

liquid  metal  condensation  interfacial  thermal 

TE 

thermoelectric 

resistance  (K/W) 

ZT 

thermoelectric  figure  of  merit 

Rte 

thermoelectric  element  radial  thermal  resistance 
(K/W) 

mean  or  average  value 

R 

gas  constant  (J/kg  K) 

Subscripts 

R 

universal  gas  constant  (J/kmol  K) 

a 

adiabatic 

s 

pool  boiling  constant  (-) 

bb 

black  body 

t 

thermosyphon  wall  thickness  (m) 

c,  C 

condenser 

T 

temperature  (K/°C) 

e,  E 

evaporator 
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HSTE 

hybrid  solar  thermoelectric  system 

t 

total 

l 

liquid  medium 

te 

thermoelectric 

p,  E 

pool  boiling  at  the  evaporator 

th 

thermal 

P 

liquid  pool 

V 

vapor  medium 

ss 

selective  surface 

vl 

liquid  vapor  interface 

sat 

saturation  conditions 

w 

wall 

(PYT)  systems  (Tripanagnostopoulos  et  al.,  2002;  Dubey 
and  Tiwari,  2008;  Garg  and  Adhikari,  1999;  Shahsavar 
and  Ameri,  2010;  Chowdhury  et  al.,  2010;  Kribus  and 
Mittelman,  2008)  use  photovoltaics  (PV)  to  generate 
power,  but  are  limited  to  low  temperature  (30-80  °C)  waste 
heat  production.  Concentrating  hybrid  photovoltaic-ther¬ 
mal  (CPVT)  systems  (Gibart,  1981;  Buffet,  1982;  Zondag, 
2005;  Coventry,  2004;  Rosell  et  al.,  2005;  Tripanagnosto¬ 
poulos,  2007)  show  higher  electrical  conversion  efficiencies 
than  PVT  systems  but  are  also  limited  to  low  temperature 
operation.  In  these  systems,  there  is  a  significant  degrada¬ 
tion  in  electrical  conversion  efficiency  when  temperatures 
increase  due  to  higher  internal  carrier  recombination  in 
the  PV  (Green,  2003;  Luque  and  Marti,  1999). 

Thermoelectrics  (TE),  in  contrast,  promise  higher  elec¬ 
trical  conversion  efficiencies  at  elevated  operating  tempera¬ 
tures  because  thermal  energy  is  directly  converted  to 
electrical  energy  via  the  Seebeck  effect.  When  a  temperature 
difference  exists  across  the  TE,  power  is  produced  with  no 
moving  parts.  Solar  TE  energy  conversion  systems  (Lenoir 
et  al.,  2003;  Omer  and  Infield,  1998;  Xi  et  al.,  2007;  Deng 
and  Liu,  2009;  Rockendorf  et  al.,  1999;  Muto  and  Chen, 
2010;  Kramer  et  al.,  2011;  Chen,  2011),  where  solar  energy 
drives  the  temperature  difference  across  the  TE,  have 
significant  potential  to  produce  electrical  energy  with 


abundant  waste  heat  to  meet  building  energy  or  IPH 
demands.  Rockendorf  et  al.  (1999)  and  Li  et  al.  (2010) 
investigated  hybrid  solar  TE  water  heaters  and  numerically 
determined  that  low  radiative  losses  and  efficient  back  side 
cooling  are  needed  to  attain  electrical  conversion  efficien¬ 
cies  up  to  30%  of  the  Carnot  efficiency.  However,  Lertsatit- 
thanakom  et  al.  (2010)  experimentally  showed  that  hybrid 
solar  TE  water  heaters  with  backside  cooling  have  efficien¬ 
cies  limited  to  0.87%  due  to  electrical  pumping  require¬ 
ments.  The  results  of  these  studies  suggest  that  the 
development  of  more  efficient  cooling  methods  for  the 
TE  back  side  is  needed  to  realize  the  potential  of  hybrid 
TE  systems.  In  addition,  parametric  optimizations  could 
extend  the  utilization  of  waste  heat  from  hot  water  systems 
to  a  variety  of  applications  that  require  a  larger  range  of 
temperatures. 

In  this  work,  we  investigated  the  design  of  a  hybrid  solar 
TE  system  where  a  thermosyphon  passively  and  efficiently 
transfers  heat  to  a  secondary  (bottoming  cycle)  application 
for  low  temperature  water  heating  or  higher  temperature 
co-generation  and  IPH.  In  addition,  thermosyphons  are 
reliable  and  can  be  used  for  a  large  range  of  temperatures 
(Faghri,  1995).  Fig.  1  shows  the  hybrid  solar  thermoelectric 
(HSTE)  configuration  we  studied  in  this  work.  A  parabolic 
concentrator  focuses  light  on  the  evaporator  section  of  the 


Condenser  Section 
Thermosyphon 

Selective  Absorber. 
Vacuum 
Glass 
Thermoelectric 
Thermosyphon 


Fig.  1.  Schematic  of  the  hybrid  solar  thermoelectric  system  (HSTE).  Solar  energy  is  focused  by  a  parabolic  concentrator  on  the  evaporator  section  of  the 
evacuated  tube  absorber  (thermosyphon),  which  heats  the  TE  hot  side.  The  resulting  temperature  difference  between  the  TE  hot  and  cold  sides  produces 
electrical  power  while  heat  carried  away  (waste  heat)  by  the  thermosyphon  (adjacent  to  the  cold  side)  is  transferred  to  the  condenser  section  for  the 
bottoming  cycle. 
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evacuated  tube  absorber  (thermosyphon),  heating  the  TE 
hot  side  located  on  the  inner  tube.  The  resulting  tempera¬ 
ture  gradient  between  the  TE  hot  and  cold  sides  produces 
electrical  power.  Meanwhile,  the  thermosyphon  interfaced 
on  the  TE  cold  side  maintains  the  temperature  and  trans¬ 
fers  the  heat  (waste  heat)  to  the  condenser  section  for  co¬ 
generation,  IPH,  or  hot  water  heating. 

We  developed  an  energy-based  model  to  investigate  the 
effect  of  solar  concentration,  TE  and  thermosyphon  mate¬ 
rial  and  geometry,  thermosyphon  working  fluid,  and  bot¬ 
toming  cycle  temperatures  on  HSTE  system  performance 
(Section  2).  The  results  from  the  model  are  subsequently 
used  to  (i)  optimize  the  HSTE  efficiency  and  radial  geome¬ 
try,  (ii)  investigate  the  effect  of  different  thermal  conductiv¬ 
ity  thermosyphon  wall  materials  to  potentially  reduce  cost, 
and  (iii)  propose  potential  commercial  applications  for 
HSTE  utilization  (Section  3).  The  outcomes  of  this  work 
aid  in  predicting  the  performance  of  five  potential  solar 
applications  and  provide  guidelines  for  the  design,  optimi¬ 
zation  and  selection  of  TE  and  thermosyphon  components 
for  future  high  performance  HSTE  systems. 


2.  Model 

2.1.  Model  formulation 

A  cross-sectional  schematic  with  the  geometric  parame¬ 
ters  of  the  HSTE  and  the  equivalent  thermal  resistance 
model  (Prasher,  2003;  Tundee  et  al.,  2010)  are  shown  in 
Fig.  2a  and  b,  respectively.  A  solar  parabolic  concentrator 
(Fig.  1)  focuses  sunlight  ( Qsoiar )  on  a  selective  surface  (SS) 
with  a  low  thermal  emissivity  and  high  solar  absorbtivity. 
The  surface  also  emits  thermal  radiation  ( Qhss )  at  a  spec¬ 
trally  averaged  emissivity  (e)  due  to  its  elevated  temperature 
(Tss).  To  accurately  capture  the  absorptive  and  emissive 
properties,  NREL  data  was  used  for  four  thermally  and 
mechanically  robust  commercially  available  surface  coatings 
(Kennedy  and  Price,  2006)  (see  Appendix  A).  The  selective 
surface  is  assumed  to  be  isothermal  along  its  length,  and 
the  temperature  drop  from  the  SS  to  the  TE  hot  side  is  across 
a  thin  film  (<500  pm)  and  therefore  can  be  neglected  (Chen, 
2011).  We  also  assumed  that  in  the  evacuated  concentric  tube 
design  there  are  only  radiative  losses  from  the  selective 
surface. 

The  net  heat  absorbed  by  the  SS  is  conducted  through 
the  TE  element  with  a  radial  conduction  thermal  resistance 
Rte,  which  leads  to  a  temperature  gradient  between  the  TE 
hot  (TJ  and  cold  ( Ts  F)  side  to  produce  TE  power  (Pte). 
The  temperature  gradient  is  dependent  on  the  TE  leg  geom¬ 
etry,  material  thermal  conductivity  ( kte ),  and  figure  of  merit 
( ZT ).  Due  to  the  high  sensitivity  of  these  parameters  on 
system  performance,  temperature  dependent  properties 
were  used  for  the  TE  thermal  conductivity  (Fig.  3a)  and  fig¬ 
ure  of  merit  (Fig.  3b)  (Minnich  et  al.,  2009;  Poudel  et  al., 
2008;  Snyder  and  Toberer,  2008;  El-Sharkawy  et  al., 
1983;  Morelli  et  al„  2008). 


An  inclined  two-phase  thermosyphon  in  contact  with 
the  TE  cold  side  transfers  heat  ( Qout )  axially  to  the  bottom¬ 
ing  cycle  application  at  a  temperature  Tc.  The  thermosy¬ 
phon  achieves  efficient  spreading  via  a  working  fluid  that 
undergoes  phase-change  due  to  the  heat  supplied  to  the 
evaporator.  The  generated  vapor  axially  transports  to  the 
condenser  section.  As  heat  is  transferred  to  the  bottoming 
cycle  application,  the  vapor  condenses  and  returns  back  to 
the  evaporator  by  gravity.  The  heat  transfer  processes  in 
the  thermosyphon  are  modeled  as  thermal  resistances 
(Fig.  2b):  R\  and  f?6  are  the  radial  conduction  wall  resis¬ 
tances  of  the  evaporator  and  condenser,  respectively;  R2 
and  R3  are  the  evaporation  and  boiling  resistances,  respec¬ 
tively;  R4  is  the  saturated  vapor  resistance  from  the  vapor 
flow  pressure  drop;  R5  is  the  condensation  resistance  at 
the  condenser;  R7  is  the  axial  conduction  thermal  resis¬ 
tance;  Rs  and  Rg  are  the  evaporation  and  condensation 
liquid-vapor  interfacial  thermal  resistances,  respectively. 
For  ideal  thermosyphon  operation,  the  temperature  drop 
from  the  evaporator  (TsE)  to  the  bottoming  cycle  applica¬ 
tion  ( Tc )  should  be  small. 


Fig.  2.  (a)  Schematic  cross-section  of  the  HSTE  system  in  a  horizontal 
orientation.  Solar  energy  ( Qsoiar)  heats  the  selective  surface  located  on  the 
TE  hot  side,  which  creates  a  temperature  gradient  across  the  TE  that 
produces  electrical  power  (P,e).  The  remaining  heat  (Q„ut)  is  transferred 
axially  by  the  thermosyphon  to  a  bottoming  cycle  application  at 
temperature  Tc.  (b)  Thermal  resistance  (TR)  model  of  the  HSTE,  where 
Rte  is  the  TE  radial  TR,  R\  and  Rg  are  the  thermosyphon  wall  radial  TR, 
R2  and  i?3  are  the  boiling  and  evaporation  TRs,  respectively,  R5  is 
condensation  TR,  R4  is  the  vapor  TR,  R7  is  the  thermosyphon  wall  axial 
TR,  and  Rg  and  Rg  are  evaporation  and  condensation  interfacial  TRs, 
respectively. 
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The  TEs  selected  for  HSTE  analysis  are  bismuth  tellu- 
ride  (Bi2Te3),  lead  telluride  (PbTe),  and  silicon  germanium 
(SiGe),  which  have  a  range  of  operating  temperatures  with 
moderate  ZTs  (Fig.  3),  relatively  low  cost,  and  commercial 
availability  (Minnich  et  al.,  2009;  Poudel  et  ah,  2008). 
Based  on  the  TE  and  bottoming  cycle  application  temper¬ 
atures  (Tc),  different  combinations  of  thermosyphon  wall 
materials-working  fluid  are  considered  that  ensure  working 
fluid  compatibility  and  high  effective  thermal  conductivity. 
Conventional  water-copper  (300-550  K),  mercury-stainless 
steel  (550-875  K),  and  liquid  potassium-nickel  (885— 
1273  K)  thermosyphons  were  investigated  for  Bi2Te3 
(300-525  K),  PbTe  (525-850  K),  and  SiGe  (850-1200  K) 
TEs,  respectively.  The  PbTe  and  SiGe  HSTEs  are  suitable 
for  medium  to  high  temperature  IPhf  applications,  while 
Bi2Te3  HSTEs  can  be  used  for  low  temperature  IPH  or  res¬ 
idential  heating.  The  TE  module  contains  negligible 


spacing  between  the  TE  legs  such  that  only  conduction  heat 
transfer  through  the  TE  legs  is  considered  (Rockendorf 
et  ah,  1999;  Li  et  ah,  2010). 


2.2.  Governing  equations 

We  solved  the  energy  equation  governing  the  thermo¬ 
electric  and  thermosyphon  over  a  wide  range  of  input 
parameters  including  solar  concentration,  TE  and  ther¬ 
mosyphon  material  and  geometry,  thermosyphon  working 
fluid,  and  bottoming  cycle  temperature.  All  transport  prop¬ 
erties  vary  with  temperature  including  the  TE  thermal  con¬ 
ductivity  ( kte ),  figure  of  merit  ( ZT ),  thermosyphon  wall 
thermal  conductivity  (kw),  and  thermosyphon  fluid 
properties  {kv,kh  pv,  ph  a,  pv,  ph  Cpv,  Cph  hfg).  Fig.  4a 
and  b  show  schematics  of  the  thermosyphon  evaporator 


(a) 


Temperature  [K] 


(b) 


Temperature  [K] 


Fig.  3.  Thermoelectric  (a)  thermal  conductivity  (kte)  and  (b)  figure  of  merit  (ZT)  as  a  function  of  temperature  for  three  TE  materials  used  for  the  model: 
bismuth  telluride  (Bi2Te3),  lead  telluride  (PbTe)  and  silicon  germanium  (SiGe)  (Minnich  et  al.,  2009;  Snyder  and  Toberer,  2008;  Poudel  et  al.,  2008). 


(a) 


I9 


- ►!  §£o  It - 


Fig.  4.  Schematics  of  the  (a)  TE  side  showing  the  film  evaporation  and  pool  boiling  region  of  the  evaporator  and  (b)  condenser  side  showing  film 
condensation  in  a  vertical  orientation.  Energy  inputs  and  outputs  are  labeled. 
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and  condenser,  respectively,  which  include  the  energy 
inputs  and  outputs  in  the  system. 

At  the  TE  side  (Fig.  4a),  the  energy  gained  from  the 
solar  heat  input  ( Qsoiar )  and  lost  from  emissive  loss  ( Q,mt) 
is  balanced  by  the  generated  TE  power  ( Pte )  and  waste  heat 
(Qour)- 

y-AcsCG  -  aBEAE(T*ss  -  T4J  =  Pte  +  Qou, 

Solar  Heat  Input(Qsolar)  Radiative  I.oss(&,«„)  TE  P'1""  Waste  Heat 

(1) 


where  C  is  the  solar  concentration  ratio,  Asc  is  the  evapo¬ 
rator  cross  sectional  area  (2 rteLe),  AE  is  the  evaporator  sur¬ 
face  area  (2n rteLe),  G  is  the  average  solar  insolation 
(1000  W/m2),  a  and  e  are  the  selective  surface  spectrally 
averaged  solar  absorbtivity  and  thermal  emissivity,  respec¬ 
tively,  oB  is  the  Stefan-Boltzmann  constant,  Tss  is  the  selec¬ 
tive  surface  and  TE  hot  side  temperature,  and  Tx  is  the 
ambient  temperature  (300  K).  The  TE  power  is  defined  as 


e  =  (Qsoiar -Q,o J  '  *l,e 

T  ~TrF 

=  (QSo,ar-QIo J-— - ^ 


(2) 


where  r\te  is  the  TE  electrical  generation  efficiency,  TsE  is 
the  TE  cold  side  temperature,  and  ZT  is  the  thermoelectric 
figure  of  merit.  The  TE  material  properties  were  averaged 


Table  1 

Thermal  resistances  for  low  and  high  temperature  HSTE  models. 
Low  temperature  model  ( T  ^  550  K) 


Rte 


(3) 


*2=2 

hr.  =  0-943[»feA^]‘  f  ^  20 

h=  A  f  >  20  (5) 

A  =  [0.997  -  O.334(coj0)o1O8][^][o-254(otiS)O'1OSi 

«3  =  0  ‘  (6) 


*5  =35^  , 

hc  =  0.943  -£c$\  *  7?  <  20 

hc  =  A  f  >  20  (8) 

A  =  [0.997  -  0.334(co^)0108][^][0'254(oti8),UO‘1 

*6=S  W 

r7  =  (10) 

R«=0  (11) 


a  See  text  for  definition  of  hEy- 
b  See  text  for  definition  of  hEp. 


over  the  operating  temperature  interval  [Tss,  Ts  E\.  At  the 
thermosyphon  condenser  side  (Fig.  4b),  the  energy  trans¬ 
ferred  from  the  evaporator  ( Qout )  is  transferred  directly 
to  the  bottoming  cycle  application  at  the  condenser  tem¬ 
perature  (Tc). 


2.3.  Thermosyphon  modeling 

The  thermosyphon  is  modeled  using  a  thermal  resistance 
network,  R  —  AT/Q,  where  R  is  the  thermal  resistance,  AT 
is  the  temperature  difference  across  each  thermal  resistance, 
and  Q  is  the  heat  transfer  (Fig.  2b).  The  resistances  are 
determined  for  low  (T <  550  K)  and  high  (T >  550  K)  tem¬ 
perature  operating  regimes  (Table  1).  At  low  temperatures, 
classical  Nusselt  condensation/evaporation  film  theory 
(Faghri,  1995;  Incropera  and  DeWitt,  1985)  is  used.  How¬ 
ever,  at  high  temperatures  with  liquid  metals  as  the  work¬ 
ing  fluid,  high  vapor  velocities  result  in  large  interfacial 
shear  stresses  and  additional  interfacial  heat  transfer  resis¬ 
tances  (i?g  and  Rg)  at  the  liquid  vapor  interface.  Therefore, 
a  modified  Nusselt  analysis  with  pool  boiling  and  thin  film 
evaporation  (Sukhatme  and  Rohsenow,  1966;  Faghri, 
1995;  Da  Cunha  and  Mantelli,  2006,  2009)  more  effectively 
captures  the  phase-change  at  the  evaporator.  In  all  cases, 
the  ratio  of  condensate  film  thickness  ((5)  and  tube  radius 
(r,)  is  assumed  to  be  very  small  (<5/r,  1),  such  that  the 


High  temperature  model  (T>  550  K) 


*«  =  5s 


(13) 


(14) 


(15)a 


i?3  = 

«4  = 


(InrtLp+izrf  )hEjl 


3  2 hcxr,Lc  J 

hc  =  0.943[^:y]8f<20 

hc  =  A  !?)!] *7?  >  20 

A  =  [0.997  -  O.334(coi0)olO8][^]IO-254(c“e,OIOSl 

r  to® 

~  2nLckw 
Rj  - 

*s  =  (j^Xi)’(f 

*See  text  for  definition  calculation  of  o' 


*9 

‘See  text  for  definition  calculation  of  o' 


(16)b 

(17) 


(18) 

(19) 

(20) 
(21) 

(22) 
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tube  can  be  modeled  as  laminar  film  condensation  on  a  flat 
inclined  surface. 

2.3.1.  Low  temperature  thermosyphon  model  ( TsE  <  550  K) 

We  considered  Nusselt  film  condensation  and  evapora¬ 
tion  for  a  specified  condenser  outer  wall  temperature  ( Tc ) 
to  obtain  the  temperature  distribution  of  low  temperature 
thermosyphons,  where  the  thermal  resistances  are  given 
by  Eqs.  (3)— (12).  For  small  thermosyphon  length  to  radius 
ratios  {Ljrt  <  20),  the  thermosyphon  tube  is  modeled  as 
laminar  film  condensation  on  a  flat  inclined  surface,  where 
the  condensed  film  returns  to  the  evaporator  by  gravity  and 
evaporates  completely  by  the  end  of  the  evaporator  section 
(i.e.,  no  liquid  pool  exists  at  the  base  of  the  evaporator). 
However,  for  larger  length  to  radius  ratios  ( Lclrt  >  20),  a 
modified  heat  transfer  coefficient  correlation  is  used  that 
better  captures  the  circumferential  liquid  film  thickness 
variation  due  to  thermosyphon  inclination  (Hussein 
et  al.,  2001).  The  vapor  is  assumed  to  be  at  saturation  con¬ 
ditions,  and  shear  forces  are  negligible  (Faghri,  1995) 
resulting  in  thermal  resistances  of  evaporation  and  conden¬ 
sation  heat  transfer  described  by  Eqs.  (5)  and  (8), 
respectively. 

To  determine  the  difference  in  saturation  temperature 
from  the  evaporator  (TsatE)  to  condenser  (Tsat  C),  the  Cla- 
peyron  relation  is  used  which  accounts  for  the  vapor  flow 
thermal  resistance  (Eq.  (4)).  While  previous  works  have 
found  this  resistance  to  be  negligible  (Prasher,  2003),  the 
moderate  heat  fluxes  and  temperatures  in  this  analysis 
can  lead  to  appreciable  temperature  drops. 

2.3.2.  High  temperature  thermosyphon  model 
(TSiE  >  550  K) 

At  higher  temperatures,  the  interfacial  shear  stresses  can 
create  significant  error  in  the  predicted  film  thickness  pro¬ 
file  using  classical  Nusselt  theory.  Therefore,  a  separate 
high  temperature  model  is  used  to  predict  thermosyphon 
performance  using  a  modified  Nusselt  model  with  pool 
boiling  heat  transfer,  where  the  thermal  resistances  are 
given  by  Eqs.  (13)-(22).  For  simplicity,  we  also  assume  that 
the  evaporator  and  wall  temperatures  are  uniform  and  use 
the  modified  Nusselt  condensation  model  in  the  condenser 
section  (Da  Cunha  and  Mantelli,  2006,  2009)  with  the  addi¬ 
tion  of  liquid  vapor  interfacial  resistances  [Rg  and  Rg). 

To  determine  the  liquid  vapor  interfacial  resistances 
(Eqs.  (21)  and  (22)),  we  assume  condensation  and  evapora¬ 
tion  coefficients  (<r')  of  0.1  which  is  appropriate  for  large 
engineering  systems  which  typically  are  difficult  to  main¬ 
tain  in  a  pure  environment  (Sukhatme  and  Rohsenow, 
1966).  The  film  evaporation  heat  transfer  coefficient  from 
the  top  of  the  evaporator  section  to  the  liquid  pool  is 

t  gcos6p,(p,  -  pv)hfg 

J'E  3  p,{Tsa,,E-Tn,E){Le-Lp) 

X  [(.4(Z,  -  Lp )  +  §1  -  A(LC  +  La)f  -  4]  (23) 


s  _  plfr/MW  -  A  _  4kip,(TsalE  -  7>) 

io  [  gcos0p,(p,  -  pv)hfg  \  ’  gcos0p,(p,  -  pv)hfg 

where  dLo  is  the  film  thickness  at  the  end  of  the  adiabatic 
section,  Lp  is  the  length  of  the  pool  region,  Le  is  the  evap¬ 
orator  length,  Lc  is  the  condenser  length,  and  L,  is  the  total 
length  of  the  thermosyphon  (Da  Cunha  and  Mantelli, 
2009).  Eq.  (23)  is  combined  with  Eq.  (15)  to  determine 
the  effective  thermal  resistance  of  the  thin  film  liquid  metal 
evaporation. 

To  incorporate  the  effect  of  pool  boiling  of  liquid  metals, 
the  heat  transfer  coefficient  is  determined  by 

hpjc  =  Cq01  P”  (25) 

where  Pr  =  P,/Pc,  C  =  13.7,  m  =  0.22  for  Pr  <  0.001,  and 
C  =  6.9,  m  =  0.12  for  Pr  >  0.001,  q  is  the  evaporator  heat 
flux,  Pc  is  the  critical  pressure  of  the  liquid  metal,  and  Pi 
is  the  liquid  pressure  in  contact  with  the  heated  surface 
(Shah,  1992).  This  average  heat  transfer  coefficient  is  used 
with  Eq.  (16)  to  determine  the  effective  thermal  resistance 
due  to  liquid  metal  pool  boiling  (R3). 

For  both  low  and  high  temperature  thermosyphon  mod¬ 
els,  we  assume  operation  below  limiting  thermosyphon 
operating  conditions.  To  verify  this  assumption,  the  con¬ 
tinuum,  sonic,  viscous,  entrainment,  and  boiling  limits 
(Faghri,  1995;  Dickerson,  1996)  are  calculated  and  com¬ 
pared  with  the  corresponding  operational  performance 
(see  Appendix  A).  In  addition,  non-condensable  gases  are 
assumed  to  be  in  negligible  amounts  as  to  not  affect  the 
heat  transfer  characteristics.  We  also  verified  that  the  film 
condensation  is  laminar  using  the  condensation  Reynolds 
number  (Re ,5  <  30).  The  transient  response  of  the  HSTE 
is  not  considered  due  to  the  relatively  low  thermal  mass 
and  quick  startup  time  of  thermosyphons  and  evacuated 
tube  solar  collectors  (Mather,  1982;  Faghri,  1995;  de  Win¬ 
ter,  1991). 

2.4.  Solution  algorithm 

The  HSTE  model  was  solved  iteratively  where  the  solar 
concentration  (C),  bottoming  cycle  temperature  (Tc),  TE 
and  thermosyphon  geometry  and  materials  and  thermosy¬ 
phon  working  fluid  are  input  parameters  to  obtain  the  tem¬ 
perature  distribution  ([T]),  thermoelectric  power  (Ple), 
waste  heat  (Qout),  emissive  loss  (Qioss)  and  system  efficiency. 
To  determine  the  transport  properties  ([R]o),  a  guessed  ini¬ 
tial  temperature  distribution  ([7]0)  is  used.  Once  initial 
properties  are  obtained,  the  model  iterates  to  determine  a 
new  temperature  distribution  ([T]])  by  solving  the  energy 
equation  (Eq.  (1))  with  thermal  resistances  (Eqs.  (3)- 
(22)).  The  transport  properties  ([E]iJ  are  recalculated  with 
the  new  temperature  distribution  and  used  in  the  next  iter¬ 
ation.  The  convergence  criterion  is  defined  as  when  the  dif¬ 
ference  between  successive  temperatures  for  each  point  is 
less  than  0.01°  (|7)+1  —  T\  <  0.01)  (see  Appendix  A).  The 
choice  of  materials  for  the  thermosyphon  and  TE  is  made 
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based  on  the  calculated  temperature  distribution  in  each 
successive  iteration. 


3.  Results  and  discussion 


Inclination  angles  from  the  vertical  (0)  up  to  30°  were 
examined  and  show  a  small  effect  on  the  performance  of 
the  HSTE  system  (<2%),  For  inclination  angles  up  to 
70°,  the  thermosyphon  can  still  operate  without  significant 
deterioration  (<10%)  which  agrees  well  with  previous 


3.1.  Model  results 

The  modeling  results  were  obtained  for  a  particular 
solar  collector  resembling  the  glass  tube  evacuated  design 
(Le  =  50  cm,  Lc  =  10  cm,  La  =  200cm,  r0  =  2.25cm, 
r,  =  2  cm,  rte  =  3  cm,  6  =  30°)  with  black  chrome  as  the 
selective  surface  which  is  stable  at  high  temperatures 
(300-800  K)  (see  Appendix  A).  Also,  the  emissivity  is  rela¬ 
tively  high  (0.08  <  e  <  0.3)  compared  to  other  selective  sur¬ 
faces,  allowing  for  a  conservative  estimate  of  performance. 
Fig.  5a  shows  the  bottoming  cycle  heat  transfer  ( Qou, )  as  a 
function  of  bottoming  cycle  temperature  (Tc)  and  solar 
concentration  (C).  Three  distinct  regimes  for  different 
TEs  exist  with  Tc.  When  the  TE  or  thermosyphon  temper¬ 
atures  exceed  550  K,  the  PbTe  TE  with  the  mercury-stain- 
less  steel  thermosyphon  replaces  the  Bi2Te3  TE  with  the 
water-copper  thermosyphon,  creating  a  discontinuity  in 
the  system  performance.  These  discontinuities  are  larger 
at  higher  Cs  because  a  larger  temperature  difference  exists 
at  higher  heat  fluxes.  A  similar  shift  occurs  at  TsJE  >  778  K 
to  a  SiGe  TE  with  the  liquid  potassium-nickel  thermosy¬ 
phon.  As  the  C  increases,  the  shift  occurs  at  a  lower  Tc 
because  the  thermosyphon  temperature  drop  is  greater, 
leading  to  a  higher  cold  side  TE  temperature.  As  a  result, 
Bi2Te3  HSTE  systems  at  high  Cs  have  a  very  narrow  oper¬ 
ation  window  (300  K  <  Tc  <  340  K). 

Fig.  5b  shows  the  emissive  loss  (Qioss)  as  a  function  of 
Tc,  C,  and  TE  material.  As  C  increases,  the  selective  surface 
temperature  ( Tss )  increases,  leading  to  higher  emissive  loss. 
Therefore,  the  decrease  in  Qout  (Fig.  5a)  is  more  pro¬ 
nounced  at  Tc  >  500  K  due  to  the  fourth  order  dependence 
of  emissive  losses  on  temperature.  Fig.  5c  shows  the  output 
TE  power  as  a  function  of  Tc,  C,  and  TE  material.  The  TE 
efficiency  is  dependent  on  the  temperature  difference 
(Tss  -  Tw  E)  across  the  TE  module  and  figure  of  merit 
(Z7)  (Section  1).  As  Tc  increases,  Pte  decreases  due  to  a 
decrease  in  Qout.  With  higher  C,  however,  a  higher  temper¬ 
ature  difference  across  the  TE  element  can  be  attained  due 
to  the  higher  heat  flux,  leading  to  a  higher  power  output 
(■ Pte )■ 

Fig.  6  shows  the  thermosyphon  evaporator  wall  temper¬ 
ature  ( Tw  E )  as  a  function  of  condenser  temperature  ( Tc) 
and  solar  concentration  (Q.  The  difference  between  evapo¬ 
rator  and  condenser  temperature  is  larger  in  the  low  tem¬ 
perature  regime  due  to  the  higher  thermosyphon  heat 
transfer  and  lower  emissive  losses.  Additionally,  low  tem¬ 
perature  copper  water  thermosyphons  have  lower  effective 
heat  transfer  coefficients  than  high  temperature  liquid 
metal  thermosyphons  due  to  the  lower  thermal  conductiv¬ 
ity,  surface  tension  and  latent  heat  of  vaporization  of 
water. 


Condenser  Temperature,  Tc  [K] 
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Fig.  5.  (a)  Thermosyphon  heat  transfer,  (b)  emissive  losses  and  (c)  TE 
power  of  the  HSTE  system  for  varying  solar  concentrations  and  condenser 
temperatures.  As  the  condenser  temperature  ( Tc )  increases,  emissive  losses 
(Qioss)  increase  while  TE  power  (P,e)  and  thermosyphon  waste  heat  (Qout) 
decrease  due  to  the  increased  surface  temperatures  (7’ss).  As  solar 
concentration  (C)  increases,  emissive  losses,  TE  power  and  heat  output 
increase. 
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experiments  and  theory  (Negishi  and  Sawada,  1983;  Hahne 
and  Gross,  1981).  Once  the  thermosyphons  are  inclined  to 
greater  than  70°,  they  cannot  be  used  for  HSTE  system 
integration  due  to  poor  performance.  Similarly,  the  ther¬ 
mosyphon  adiabatic  section  length  ( La )  shows  minimal 
effect  (<0.1%)  on  performance  at  low  condenser  tempera¬ 
tures  (Tc<  550  K)  and  solar  concentrations  (C  <  50),  indi¬ 
cating  that  the  saturation  temperature  drop  associated  with 
the  vapor  pressure  drop  in  the  thermosyphon  is  negligible 
in  these  regimes.  However,  as  discussed  in  Section  2.3, 
the  saturation  temperature  drop  needs  to  be  considered 
in  liquid  metal  HSTE  systems  due  to  a  significant  vapor 
pressure  drop  at  high  solar  concentration  ratios. 

3.2.  HSTE  optimization 

To  appropriately  select  the  highest  performance  solution 
for  a  given  application,  optimization  of  the  HSTE  system  is 
needed.  The  HSTE  is  optimized  based  on  a  combination  of 
two  efficiencies:  (1)  the  electrical  efficiency  (r],e)  defined  as 
the  TE  efficiency  (Eq.  (2))  and  (2)  the  thermal  efficiency 
( rfc )  from  waste  heat  ( <2/avs)  defined  using  an  exergetic 
approach  with  the  Carnot  efficiency  (Eq.  (26))  (Rosen 
et  al.,  2004) 


The  defined  thermal  efficiency  serves  as  an  ideal  upper 
limit  to  the  amount  of  work  obtained  from  the  waste  heat 
if  a  heat  engine  is  used  as  the  bottoming  cycle.  The  overall 
system  efficiency  {tihste)  is  therefore  defined  as  the  ratio  of 
useful  energy  extracted,  including  the  thermal  {Qout)  and 
electrical  components  ( Pte ),  compared  to  the  total  energy 
input  (see  Appendix  A  for  full  derivation) 
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Fig.  6.  Thermosyphon  evaporator  temperature  for  varying  solar  concen¬ 
trations  and  condenser  temperatures.  As  the  condenser  temperature  (Tc) 
increases,  emissive  losses  (2/„sJ  increase  while  thermosyphon  waste  heat 
(Qout)  decreases  resulting  in  a  smaller  temperature  drop  across  the 
thermosyphon.  As  solar  concentration  (Q  increases,  the  heat  output 
increases,  resulting  in  a  larger  temperature  drop.  Additionally,  liquid 
metal  thermosyphons  show  smaller  temperature  drops  due  to  more 
efficient  heat  transfer. 


_  Useful  Energy  Out 
n"sn:  =  Incident  Solar  Energy 

=  (>7te  +  »lc-Mc)(l-|f)  (27) 

Fig.  7  shows  the  HSTE  efficiency  as  a  function  of  bot¬ 
toming  cycle  temperature  and  solar  concentration.  Similar 
to  Fig.  5,  discontinuities  due  to  temperature  operation  lim¬ 
its  of  the  TE  and  thermosyphon  exist  (see  Section  3.1).  To 
obtain  HSTE  efficiency  for  a  particular  bottoming  cycle 
application,  the  end  use  temperature  at  which  the  thermal 
energy  will  be  transferred  to  must  be  known.  For  example, 
if  the  HSTE  system  is  to  be  used  for  space  heating,  the  tem¬ 
perature  is  set  to  the  condenser  temperature  ( Tc ),  and  the 
system  performance  can  be  obtained  from  Fig.  7. 

3.2.1.  Effect  of  bottoming  cycle  temperature  (Tc)  and 
concentration  (C) 

Fig.  7  shows  that  the  HSTE  efficiency  {y\hste)  has  opti¬ 
mal  values  as  a  function  of  bottoming  cycle  temperature 
(T<:).  The  initial  increase  is  due  to  the  increase  in  thermal 
efficiency  with  increasing  temperature.  However,  as  Tc  con¬ 
tinues  to  increase,  the  surface  temperature  (Tss)  of  the  TE 
element  reaches  a  point  where  the  emissive  losses  ( Qioss ) 
begin  to  dominate.  As  a  result,  with  any  additional  increase 
in  Tc,  the  efficiency  decreases  due  to  the  fourth  order  tem¬ 
perature  dependence  of  the  emissive  losses.  At  constant  Tc, 
the  system  efficiency  ( Ihste )  increases  with  increasing  solar 
concentration  (Q  due  to  higher  heat  transfer  (Qout) 
through  the  thermosyphon  and  greater  thermal  efficiency. 
Concentrations  beyond  100  suns  may  be  more  advanta¬ 
geous  but  may  have  economic  implications  on  the  con¬ 
struction  of  the  concentrator;  and  therefore  were  not 
considered  here. 


Condenser  Temperature,  Tc  [K] 


Fig.  7.  Efficiency  of  the  HSTE  system  for  varying  solar  concentrations  (C) 
and  bottoming  cycle  temperatures  (Tc).  Optimal  system  efficiencies  exist 
which  balance  the  thermal  efficiency  and  emissive  power.  Increasing  the 
solar  concentration  also  increases  efficiency  due  to  a  higher  energy  input 
and  thermal  efficiency. 
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3.2.2.  Effect  of  TE  leg  length  (Lte)  and  thermosyphon  size 
(rjrte) 

Fig.  8a  examines  the  HSTE  system  efficiency  with  TE  leg 
length  (Lte)  and  cross-sectional  radii  ratio  ( rjrte )  for  a  par¬ 
ticular  solar  concentration  (C=50)  and  bottoming  cycle 
temperature  ( Tc  =  700  K)  determined  from  Fig.  7.  As  Lu. 
increases,  the  system  efficiency  decreases  due  to  the  addi¬ 
tional  thermal  resistance  of  the  TE  leg,  leading  to  an  ele¬ 
vated  surface  temperature  ( Tss)  and  higher  emissive  loss 
( Qioss )•  The  TE  leg  length,  however,  has  different  effects 
on  system  performance  depending  on  the  bottoming  cycle 
temperature  ( Tc ).  At  high  Tc,  small  increases  in  Lte  result 
in  larger  decreases  in  system  efficiency  due  to  the  higher 
emissive  losses  at  increasing  temperatures.  However, 
Fig.  8b  shows  that  the  TE  power  increases  with  increasing 
TE  leg  length  because  collector  area  increases  and  a  higher 
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Fig.  8.  Optimization  results  for  a  PbTe  HSTE  at  C  =  50  and  Tc  =  700  K 
showing  (a)  system  efficiency  and  (b)  TE  power.  An  increase  in  TE  leg 
length  (Lte)  decreases  efficiency  and  increases  TE  power  due  to  a  larger  TE 
thermal  gradient.  As  the  TE  leg  length  increases,  the  maximum  operating 
temperature  (776  K)  is  reached,  at  which  point  the  performance  decreases 
to  zero  (gray  area).  In  addition,  for  small  thermosyphon  radii  (bottom 
white  area),  heat  pipe  limitations  (e.g.,  sonic  limit)  prohibit  operation. 


TE  temperature  gradient  exists.  Depending  on  the  power 
needs  of  the  application,  increasing  TE  leg  length  may  have 
advantages.  For  example,  if  the  application  has  larger  elec¬ 
trical  demands,  it  may  be  more  favorable  to  sacrifice  over¬ 
all  system  efficiency  for  electrical  production. 

Fig.  8  also  shows  that  as  the  radial  ratio  of  the  ther¬ 
mosyphon  ( rjrte )  decreases  for  a  constant  LTE,  the  system 
efficiency  decreases  due  to  reduced  area  for  heat  transfer 
through  the  thermosyphon  and  the  TE  power  decreases 
due  to  reduced  selective  surface  area  for  solar  input.  As 
the  TE  leg  length  increases,  the  maximum  operating  tem¬ 
perature  (776  K)  is  reached,  and  the  performance  decreases 
to  zero  (gray  area).  In  addition,  for  small  thermosyphon 
radii  (bottom  white  area),  heat  pipe  limitations  (e.g.,  sonic 
limit)  prohibit  operation.  Furthermore,  as  the  thermosy¬ 
phon  radius  becomes  less  than  r„  ~  2  mm,  the  model  is 
no  longer  accurate  because  the  thermosyphon  pipe  wall 
cannot  be  modeled  as  a  flat  plate.  Similar  investigations 
on  geometry  were  performed  for  other  optimal  C  and  Tc 
values,  and  show  comparable  results  (see  Appendix  A). 

3.2.3.  Effect  of  thermosyphon  material 

In  an  effort  to  broaden  thermosyphon  material  selection, 
we  investigated  the  effect  of  thermosyphon  wall  thermal 
conductivity  ( kw )  on  performance.  Only  the  low  tempera¬ 
ture  model  was  considered  due  to  material  compatibility 
issues  for  high  temperature  liquid  metals.  Fig.  9  shows 
the  HSTE  system  efficiency  ( rj hste)  as  a  function  of  solar 
concentration  (C)  and  thermosyphon  wall  thermal  conduc¬ 
tivity  (kw)  for  a  Bi2Te3  TE  at  Tc  =  470  K.  As  kw  decreases 
from  10  W/mK,  the  efficiency  does  not  decrease  apprecia¬ 
bly,  which  indicates  that  the  radial  conduction  resistance 
of  the  thermosyphon  wall  is  not  dominant.  However,  as 
the  thermal  conductivity  decreases  below  approximately 
1.2  W/mK,  system  efficiency  begins  to  decrease  more 


kw  [W/mK] 

Fig.  9.  HSTE  system  efficiency  for  varying  solar  concentrations  (C)  and 
thermosyphon  wall  thermal  conductivities  (kf)  for  a  Bi2Te3  TE  at 
Tc  =  470  K.  if  is  the  HSTE  efficiency  at  high  thermal  conductivities 
(if  =  UhsteI^w  =  10  W/mK)),  which  asymptotes  to  a  constant  value.  For 
solar  concentrations  below  4  suns,  materials  with  thermal  conductivities 
larger  than  1 .2  W/mK  have  comparable  system  performance 
(if  ~  Uhste)- 
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Table  2 

HSTE  potential  applications  and  performance  (per  unit  length  of  evaporator  section). 


HSTE  application 

^7 HSTE  (%) 

77  (K) 

C 

Pte  (W/m) 

Qout  (kW/m) 

Residential  heating 

15.2 

360 

50 

200 

4 

Solar  air  conditioning11 

25.4 

400 

50 

100 

4 

Low  temperature  IPH  -  chemical  dryingb 

34.4 

500 

100 

250 

7.6 

Medium  temperature  IPHC 

48.1 

700 

100 

140 

7 

High  temperature  IPH  -  aluminum  smeltingd 

52.6 

776 

100 

60 

7.5 

a  Clausse  et  al.  (2008),  Desideri  et  al.  (2009),  El  Fadar  et  al.  (2009),  Fong  et  al.  (2010),  and  Lecuona  et  al.  (2009). 
b  Demeter  et  al.  (1991)  and  Kalogirou  (2003). 
c  Grass  et  al.  (2004)  and  Kulkarni  et  al.  (2008). 
d  Murray  (1999). 


significantly  from  the  efficiency  at  high  kw  values  (rj*).  The 
results  indicate  that  materials  such  as  glass  can  be  used  for 
thermosyphons  in  HSTE  systems  when  solar  concentra¬ 
tions  are  below  4  suns,  which  can  reduce  material  costs 
(Window  and  Harding,  1984).  A  glass  thermosyphon 
design  can  be  integrated  into  a  bottoming  cycle  process 
via  a  header  and  manifold  design  most  commonly  seen  in 
current  solar  thermal  applications  (Zhiqiang,  2005;  Win¬ 
dow  and  Harding,  1984). 

To  achieve  the  optimal  HSTE  system  for  a  prescribed 
bottoming  cycle  using  the  framework  described  above, 
the  following  procedure  should  be  followed:  (1)  An  initial 
geometry  should  be  specified  for  which  an  optimum  solar 
concentration  and  bottoming  cycle  temperature  can  be 
determined  (e.g.,  Fig.  7).  (2)  The  thermosyphon  and  TE 
geometry  should  be  optimized  based  on  results  from  step 
1  (e.g.,  Fig.  8).  (3)  A  new  geometry  can  then  be  selected 
that  meets  the  power  and  heat  requirements  of  the  partic¬ 
ular  application.  (4)  Steps  1-3  need  to  be  repeated  until 
convergence  is  reached.  In  practice,  system  constraints  such 
as  the  bottoming  cycle  temperature  and  solar  concentra¬ 
tion  are  specified,  which  simplifies  the  optimization  routine 
to  geometry  only.  Step  4  was  not  needed  for  this  study 
because  the  geometry  of  interest  was  specified. 

Five  applications  for  HSTE  systems  are  proposed  in 
Table  2  showing  applications  requiring  relatively  high  heat 
output  compared  to  electrical  power  production,  such  as 
residential  heating,  solar  AC,  and  industrial  process  heat¬ 
ing,  are  ideal  for  HSTE  integration.  The  results  also  dem¬ 
onstrate  the  applicability  of  HSTEs  for  high  temperature 
applications  such  as  chemical  drying  and  aluminum  smelt¬ 
ing,  previously  not  possible  with  existing  PVT  systems. 

4.  Conclusions 

A  hybrid  solar  thermoelectric  system  using  a  thermosy¬ 
phon  to  efficiently  transport  heat  for  a  bottoming  cycle 
(HSTE)  was  studied  over  a  wide  range  of  temperatures 
(300-1200  K),  solar  concentrations  (1-100  suns),  as  well 
as  thermosyphon  and  TE  materials  and  geometries.  An 
energy-based  model  for  the  TE  and  thermosyphon  was 
developed  to  predict  the  temperature  distribution  to  deter¬ 
mine  the  performance  of  the  HSTE  system.  The  results 
show  that  HSTE  system  efficiency  is  a  strong  function  of 


solar  concentration  and  bottoming  cycle  temperature  due 
to  the  coupling  between  temperature  and  thermosyphon 
performance,  TE  performance  and  selective  surface  proper¬ 
ties.  As  the  bottoming  cycle  temperature  increases,  the 
thermal  efficiency  increases  to  an  optimum  critical  temper¬ 
ature.  Beyond  this  temperature,  the  emissive  losses  domi¬ 
nate,  resulting  in  a  decrease  in  HSTE  efficiency.  As  solar 
concentration  increases,  both  thermal  and  thermoelectric 
efficiencies  increase  due  to  higher  heat  fluxes  to  the  HSTE. 
Geometric  optimization  of  the  HSTE  also  shows  at  higher 
TE  leg  lengths  and  radii  ratios,  higher  levels  of  electrical 
power  and  waste  heat  could  be  obtained  but  with  decreased 
efficiency.  A  range  of  optimum  efficiencies  were  deter¬ 
mined,  the  highest  of  which  include:  34.4%  (Tc  =  500  K, 
C  =  50),  48.1%  (Tc  =  700  K,  C=100),  and  52.6% 
(Tc  =  776  K,  C  =  100).  The  results  from  varying  the  ther¬ 
mosyphon  wall  material  show  when  wall  thermal  conduc¬ 
tivities  exceed  1 .2  W/mK,  system  efficiency  is 
approximately  constant,  indicating  glass  thermosyphons 
could  be  used  for  low  temperature  ( T  <  550  K)  HSTEs  to 
potentially  reduce  material  cost.  The  outcomes  of  this 
study  indicate  a  major  benefit  of  using  the  HSTE  system 
where  unlike  PVT  systems,  they  can  be  utilized  at  high  tem¬ 
peratures  while  maintaining  electrical  conversion  efficien¬ 
cies  comparable  to  or  greater  than  room  temperature 
operation.  The  development  of  the  HSTE  system  extends 
the  applicability  of  hybrid  solar  thermal  to  higher  temper¬ 
ature  processes  such  as  chemical  drying  and  aluminum 
smelting,  which  mainly  require  high  quality  heat.  This 
study  serves  as  a  framework  for  selection  and  optimization 
of  HSTE  system  configuration  based  on  the  end  use 
application. 
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